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ABSTRACT 



It is well known that processes of first order phase transitions may have occurred in the 
inflationary era. If one or more occurred well before the end of inflation, the nucleated bubbles 
are stretched to large scales and the primordial power spectrum contains a scale dependent 
non-Gaussian component provided by the remnants of the bubbles. We predict the anisotropics 
in the cosmic microwave background (CMB) induced by inflationary bubbles. We build a general 
analytic model for describing a bubbly perturbation; we evolve each Fourier mode using the 
linear theory of perturbations from reheating until decoupling; we get the CMB anisotropics by 
considering the bubbly perturbation intersecting the last scattering surface. The CMB image of 
an inflationary bubble is a series of concentric isothermal rings of different color (sign of ST/T) 
on the scale of the sound horizon at decoupling (< 1° in the sky); the resulting anisotropy is 
therefore strongly non-Gaussian. The mean amplitude of ST/T for a bubble of size L follows 
the known estimates for linear perturbations, ST/T ~ Sp/ p ■ (L / H^ 1 ) 2 . In particular, bubbles 
with size corresponding to the seeds of the observed large scale voids (tens of comoving Mpc) 
induce an interesting pattern of CMB anisotropics on the sub-degree angular scale, to be further 
investigated and compared with the forthcoming high resolution CMB maps provided by the 
MAP and the Planck experiments. 



Subject headings: 
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Introduction 



Inflation is a phase of accelerated expansion driven by very high energy quantum fields (kT ~ 10 15 
GeV), which leaves observable cosmological traces at the present. The perturbations we observe today in 
the universe are therefore the most direct link we have to very high energy physics. Today, this fundamental 
topic is under investigation from several points of view: the perturbations imprinted on the cosmic 
microwave background (CMB) at the recombination era, the possible traces on the cosmic background 
of gravitational waves, the direct reconstruction of the three dimensional matter distribution traced by 
galaxies and their peculiar velocities in the modern redshift surveys and others. 



The existence of large voids in the galaxy distribution, known for a long time (Vogcley et al. 1991 
Kauffmann & Fairall, 1991), is now being established quite firmely. Recently, several almost sphe rical 

l 



Mpc have been detected and catalogued in the SSRS2 survey ( El- Ad, 

lOh- 1 



voids with mean diameter of 40/i 

Piran & da Costa, 1996); the surrounding walls have a thickness smaller than 10/i -1 Mpc; the largest 
voids reach the spatial limits of the survey and their fraction of volume is about 40%; if smaller voids are 
also considered the voids volume fraction has a lower limi t of 75%. These amazing results are confirmed 
by the analysis on IRAS (El- Ad, Piran & da Costa, 1997); these papers contain very impressive pictures 
reproducing the large voids in a three dimensional view. Although rare isolated fai nt galaxies arc found 
inside the voids, preliminary analysis of the peculiar velocities of the walls galaxies (Da Costa ct al. 1996) 
indicate the almost complete absence of matter (baryonic and dark) in the central cavities, t he matter 
is therefore confined between the voids, and this su rely affects the cluster distribution, that recently has 



been found more correlated than expected (see e.g 

matter distribution on ~ lOO/i" 1 Mpc is revealing an evidence of extra power (Landy et al 



Einasto J. et al. 1997 and references ther ein); also, the 

1996|). The 



observed inhom qgeneities constrain t he minimal inflationary scenario (slow rolling and purely Gaussian 



fluctuations, see Padmanabhan 199 



under investigation (Rctzlatt et al 



and the consequences upon the viable parameter space are currently 
[997 | and references therein). Topo logical defects provided mainly by 



posf-inflationary phase transitions (Pen ct al. 1997, Durrer et al. 1996 and references therein), have been 
proposed to be the seeds for the observ ed inhomogcncitics. and their CMB imprints have been predicted 
(Turok et al. 1997, Durrer et al. 1996); h owever, recent results revealed some power deficit both on the 
CMB and on the perturbations spectrum (Albrecht et al. 1997). A distinctive feature of these models 
is a non-Gaussian perturbation spectrum; non-Gaussianity means that the phases of the components of 
the perturbations are correlated in some way. Getting ready to detect eventual non-Gaussianity in the 
CMB high resolution maps pr ovided by the next experiments is a non-tr ivial theoretical challenge that is 
currently under development (Ferreira & Magueijo 1997, Amendola 1996 and references therein). 



This paper aims at predicting the CMB anisotropies generated by models of first order inflation in 
which the perturbation spectrum is in general the sum of the Gaussian slow rolling fluctuations and the 
scale dependent non-Gaussian perturbations from the remnants of nucleated bub bles. First order phase 
transitions in the early universe are predicted by many inflationary theories (see [Kolb 1991 ), and one of 
the most interesting ideas introduced in inflationary cosmology in recent years is tFiat in general the phase 
transition occurs during the inflationary slow rolling (Occhionero & Amendola 1994, Amendola et al. 1996). 
Like in ordinary first order inflation, the nucleated bubbles grow as the causal horizon until the collision 
with other bubbles stops the overcomoving growth; however in these models the true vacuum energy is not 
zero, thus the bubbles may be in the linear regime and, most important, the amount of e— foldings between 
the phase transition and the end of inflation allows the bubbles to grow comovingly. This phenomenology 
has very interesting consequences for what concerns the gravitational radiation produced by the colliding 
bubble s: the characteristic pea k frequency, exponentially redshifted by inflation, may fall in a detectable 
range (Baccigalupi et al. 1997). At reheating, the energy stored in the bubble wall is converted in matter 
and radiation, leaving a bubbly trace in the cosmic density. In the case in which a phase transition occurred 
at about 50 e-foldings before the end of inflation, the remnants of the bubbles have comoving sizes of tens 



of Mpc; they could be the seeds of the large voids observed today, a possibility firstly suggested in La 1991 
and then further investigated (see e.g. Amendola et al. 1996, Occhionero et al. 1997). 



In a previous work ( Baccigalupi, Amendola fc Occhionero 1997| ) we computed the CMB anisotropies 
induced by non-linear voids and we obtained from COBE data an upper limit of about 100/t -1 Mpc to t he 



present radius of a non- linear void sitting on the last scattering surface (LSS); in Occhionero et al. 1997 



considered mainly the time evolution of spherical bubbles in the matter dominated era, as the outcome of 
linear ones at decoupling and also we obtained (again from COBE) an upper limit of about 100/i -1 Mpc to 
the present radius of linear bubbles on the LSS. Here we compute the detailed CMB anisotropies generated 
by an inflationary bubble. We build an analytic model for describing a general linear bubbly perturbation; 
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we evolve each Fourie r mode from reheating to decoupling using the theory of linear perturbations (see 
Hu & Sugiyama 1995 and references therein). Finally we get the CMB anisotropics by considering the 



intersection between the bubble and the LSS. We shall focus the analysis on bubble comoving sizes of tens 
of Mpc at decouplin g (angular scale below 1°). since they correspond to the observed structures and voids 
recently catalogued ([El-Ad, Piran fc da Costa, 19961, 1997). 



In Amcndola et al. 1997 we consider a realistic CMB fluctuation field produced by a primordial power 
spectrum containing ordinary Gaussian fluctuations and a non-Gaussian bubbly component, as predicted 
by the first order inflation models mentioned above; the volume fraction (~ 50% at the present) occupied 
by the bubbles allow them to significatively affect the large scale structure: the impact of this model on the 
CMB angular power spectrum is predicted. 

Our computations shall be compared with the high resolution CMB maps provided by the next MAP 
and Planck experiments, that will provide the whole sky CMB anisotropy maps down to 5T/T ~ 10 -6 and 
with angular resolution of 10 . 

The paper is organized as follows: section II is dedicated to the analytical construction of the 
perturbation field; in Section III we apply the linear perturbation theory and evolve the bubbly perturbation 
from reheating until decoupling focusing on the corresponding perturbation in the photon baryon plasma; 
in Section IV we compute the bubbly CMB anisotropics; finally, Section V contains the conclusions. 



2. A general analytic model 



In principle, it is possible to calculate the exact bubble shape by considering a specific first order 
inflationary model and the fields involved in the process. However, since a bubbly density perturbation is 
specified by few parameters, we choose here to parametrize the inflationary bubble directely at reheating. 
Bubbles from a phase tran sition occurred early during inflation are stretched co movinglv to large scales 
by the remaining inflation (Occhionero & Amendola 1994, |Amendola et al. 199"6 ). At reheating they leave 
a perturbation in the cosmic density. A compensated spherical bubble is described by its central density 
contrast 5, comoving radius R and thickness a of the compensating shell: 



5 = 



R , a = 



AR sh 



ell 



r=0 
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(1) 



Defining a radial coordinate x = r/R in the comoving gauge, we model naturally two zones in a bubble, the 
central underdensity and the compensating shell; the width of the latter is divided in inner (<7j) and outer 
(cr ), with <7j + <j = a. We model the bubble as follows: 



Sp 
P 



-5 (x<\ — (Ti) , 



— = A + B sin 
P 

5 -l = C + Dsin 
P 



7T / Oi 

a \ 2 



(1 - (Ti < X < 1) , 

(1 < x < 1 + a ) . 



(2) 



The constants are determined by the requests of continuity and compensation: 
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The Fourier transform of the above density contrast in a volume V is 
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This model allows for the description of any compensated spherical underdensity with the choice of 
parameters 5, R and ai,a . Figure (1) shows the radial density contrast (top) and the corresponding Fourier 
transform (bottom) for different shapes. Note that the latter is strongly scale dependent: it is characterized 
by a peak at k ~ 2tt/2R, corresponding to the diameter of the bubble. On smaller wave numbers, the 
spectrum falls as k, while the oscillations at higher k describe the compensating shell. Left panels display 
the dependence on the shell width by maintaining cr^ — cr ; right panels display the analysis of asymmetric 
shells <Ji ^ a . The variation of width and shape of the compensating shell does not affect substantially 
the amplitude and the oscillating behaviour of the Fourier transform. This makes the results of the next 
Sections almost independent of the exact shell shape; for this reason we put the emphasis on the more 
interesting physical observables R and (5, and simply assume a = 2<r, = 2a D = .3 (the solid lines in Fig.(l)). 
Also, independent reasons to neglect the dependence on the shell's shape arise from CMB damping of 
anisotropies on very small angular scales (Silk damping), to be described in the next Section. Finally, note 
that the perturbation (and its CMB counterpart) scales linearly with 8, that has been factored out in the 
axis labels. 

In the next Section we shall consider the perturbation described here in the cosmic medium at reheating 
and we will evolve it until decoupling. 



3. From reheating to decoupling 

Here we want to compute in detail the CMB perturbations induced by the remnants of inflationary 
bubbles from reheating until decoupling, leaving to the next Section the analysis of the corresponding 
CMB anisotropies. We p erform the computations in a standard scen ario (flat CDM, low baryon density 



(flh — .06). h—.5). As in Baccigalupi, Amcndola fc Occhionero 1997 we adopt the approach developed in 
Hu fc Sugiyama 1995|, since it is applicable to any kind ot scalar linear perturbation and very accurate 



(5 T rD%J! We report here only the basic and intuitive features of the method. 

Let us consider the Fourier modes (Q) as initial condition in the matter distribution at reheating. 
Let the scale factor be a(rf) and the time variable be the conformal time r\ — L dr/a(r),. Each mode 
(8p/ ' p){k,rj) evolves according to the linear theory of perturbation and its gravitational potential induces 
the corresponding time evolving perturbation (ST/T)[k,r)) in the photon baryon plasma. At decoupling, 
the latter imprints anisotropies (8T/T)(9,<f>) in the CMB. Summarizing, the process involves the following 
steps: 

So St St St 

— (M)^— (M) revolution -> —(k,T} D ) ^ — (0,0) • (5) 
pi 11 

The photon-baryon plasma contains zero-point adiabatic (ST/T)o(k,rj) and pure velocity V(k, rj) 
perturbations; at decoupling the Sachs- Wolfe effect (SW, (8T/T)sw{k, J7)) given simply by the gravitational 
potential, must be added. In the tight coupling approximation, the equations driving the time evolution are 
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(Bu fc Sugiyama 1995Q 
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the little difference between $ and 'J', due to the radiation components of the cosmic fluids, is taken into 
account in the computations. Below the Silk damping scale, the photon baryon plasma perturbations are 
damped by a factor exp[— k 2 /k^ (?])], where the damping scale fc^, 1 is of the order of ten comoving Mpc at 
decoupling and soon there after grows to the horizon scale. We simply put the spectrum (ji|) as the initial 
density perturbation; through linear theory, encoded for the CMB in (pR), we follow its evolution until 
decoupling. 

The above procedure allows to perform the first three steps indicated in (^); however, since the 
perturbations here are non-Gaussian, the computation of the CMB angular anisotropies can not be done 
automatically as for Gaussian perturbations, but requires some more details, that we leave to the next 
Section. Note that the above procedure is general and does not depend on the particular initial perturbation 
spectrum; therefore the CMB perturbation resulting from the evolution of any non-Gaussian linear feature 
from reheating may be calculated from its Fourier transform using the above procedure. Also, because 
of linearity, the time evolution of our bubbles or any kind of non-Gaussian perturbations present in the 
cosmic fluid at reheating does not depend on the presence of another (e.g. Gaussian) component in the 
perturbation spectrum. This is of particular interest for first order inflation models which provide two 
independent perturbation mechanisms, (bubbles and Gaussian fluctuations from the slow rolling field in the 
present case) as already mentioned in the Introduction. 

Let us concentrate now on the CMB perturbations induced by bubbles. We evolve the photon baryon 
perturbations until decoupling, when we antitransform to get the relevants effects as a function of r. Figure 
(2) shows the adiabatic (top), Doppler (middle) and SW effect (bottom). The bubbles under investigation 
are charecterized by a central density contrast 8 at decoupling that we put on the y-axis units, and by 
R = 10, 15, 2 Qh~ 1 Mpc; these radii are compa rable with the size of the observed large voids in the galaxy 
distribution (El-Ad, Piran & da Costa, 1996,1997). The most interesting feature of the results shown in 
Fig. (2) can be noted by looking at the curves and the scales: the SW effect strictly does not extend beyond 
the size of the perturbation, while the adiabatic and Doppler ones are not vanishing until 60 -j- 80/i _1 Mpc 
from the center of the bubble, independently of its radius. This is infact the amplitude of the sound 
horizon of the photon-baryon plasma at decoupling. Just like a pebble in a pond, the initial small bubbly 
perturbation is propagating beyond the initial scale, reaching the scale of the sound horizon at the time in 
which we are examinating it. A sound wave with negative and positive crests is well evident at the sound 
horizon position, both in the adiabatic and Doppler effects. We have further investigated the motion of 
this wave for the case R — 20/i _1 Mpc. Fig. (3) shows in detail the wave at different times. Two physically 
expected effects are immediately evident: the wave is travelling outward the bubble and simultaneously 
damped due to the Silk damping, that is becaming more and more active with time. 



Also note as the SW effect (and roughly the adiabatic one) follows the known amplitude ( Padmanabhan 
1993) for linear perturbations ST jT ~ —S(R/H^ 1 ) 2 . We have not found a significative dependence on the 

shell's shape: for instance, by taking a = a/2 and repeating the same computations, the results change of a 
few percent; this is due to the Silk damping, the effect of which is to erase the small scale CMB anisotropies 
at decoupling; in fact, the damping exponential mentioned above is 1/e at about IO/1 -1 Mpc in the present 
CDM model; thus, the damping is particularly effective on the small scales involved by the shell (a few 
comoving Mpc). 
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As we will see in the next Section, the above results imply that the CMB anisotropics of an inflationary 
linear bubble consist in general in a few concentric rings of alternate colour (sign of ST/T) on the sub-degree 
angular scale. 



4. Anisotropies 

The perturbations in the photon-baryon plasma are photographed at last scattering and described by 
the CMB anisotropy field: 

^(M)=J>JmlWM) ■ (9) 
Im 

Concerning the anisotropy amplitude, the physical information is contained in the correlation function 

W 1 / n V 1 / n' / nn'=cos0 J 47F 

where Pi are the Legendre polynomials and the average is taken over all the sky directions n and n' 
separated by the angle 9. The coefficients J2m \ a im\ 2 are conventionally renamed as (21 + 1)C/ in a Gaussian 
theory; for the functional properties of the P i, the I— th coefficient describes the CMB anisotropies at the 



angular scale 9 — n/l (see White et al., 1994 and references therein). The coefficients Ci of the expansion 
in Legendre polynomials do not determine uniquely the anisotropies, but only their amplitude at a given 
angular scale specified by I; they completely neglect the phases of the ai m coefficients, and therefore describe 
the anisotropies only if the phases are randomly distributed (Gaussian perturbations). If on the contrary 
the anisotropies are the outcome of some coherent, i.e. non-Gaussian, structure, th e Ci are indicative of 
the amplitude only, but not of the coherence. This point is of particular interest in [Amcndola et al. 1997 . 
where the subject is the impact of a bubbly distribution on the Ci spectrum. 

Here we want to compute the CMB anisotropies (|9|) generated by remnants of inflationary bubbles. On 
a particular direction n, the anisotropy is given by 

where P(rj) is the probability for a photon to be last scattered between times rj and r\ + drj; (ST /T) n (rj) 
is the CMB perturbation field at causal distance i]q ~ rj on the direction n. The perturbations examined 
here are spatially localized: they have a center and vanish beyond a sound horizon distance from it. 
Consequentely, in performing the computation (pT|), we must specify the position of the bubble's center 
with respect to the LSS. We choose to orient the polar axis towards the center of the bubble, so that due 
to its spherical symmetry, the anisotropy field (m is a function only of the polar angle 9, that is the angle 
between the bubble center direction and the particular line of sight n; the same angular choice was made in 



Baccigalupi, Amendola fc Occhionero 1997 . The remaining degree of freedom regards the distance of the 
bubble on the line of sight; its perturbation on the CMB depends on the distance of the bubble's center 
from the LSS (hereafter d) and vanishes of course if the perturbing zone does not intersects the LSS: to 
exemplify this dependence, we shall display in the sequel three interesting cases: d = —R, and +R. Let us 
consider a photon last scattered at some rj inside the bubble and running towards the observation point on 
some direction 9; the ST/T that it carries depends on several variables: 9, d and r\ uniquely determine the 
radial coordinate (radial with respect to the bubble's center of course) at which it decouples; we compute 
the adiabatic, SW and Doppler effects and the overall CMB perturbation 

ST ^ djr))= ( 5 J.) + ( S _T) +V .cosa , (12) 

SW 



T \ T n \ T 



where a is the angle between the radial direction and the line of sight, and of course accounts for the 
direction dependence of the Doppler effect. The above expression means that for each rj we antitransform 
(ST/T)(k,rf) as obtained from the procedure of the previous Section. Note that this is particularly accurate, 
since the behaviour of the Silk damping scale (that suddenly after the decoupling grows up to the horizon 
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scale) is taken into account. Finally we get the CMB anisotropics simply averaging on all 77's, weighted 
with the last scattering probability P(rj) (|Tl|). The integrated SW effect, due to the time change of the 
gravitational potential while the ph oton crosses the bubble, is not taken into ac count since it is completely 



negligible, as in the non-linear case (Baccigalupi, Amendola & Occhionero 1997) 



Figure (4) reports the CMB anisotropics from bubbles lying on the LSS; the radii are the same as in 
Fig. (2). In each panel, three different lines show the interesting variation of the angular dependence of 
8T/T with the relative positions of the LSS with respect to the bubble's center: d = —R, 0, +R (respectively 
short dashed, solid and long dashed). Of course the whole signal scales with 8 and grossly the features 
of Fig. (2) are respected. Note the effects of the variation of d: the negative strong spot near the center, 
coming mainly from the adiabatic effect (see also Fig. (2)), is depressed if the bubble is displaced behind or 
beyond the LSS. 

Finally in Figure (5) the CMB sky contains only one bubble with R = 20h~ 1 Mpc and 8 — .005 (top), 
8 = .01 (bottom), lying on the LSS, corresponding to the solid line on the bottom panel in Fig. (4); in the 
same sky region we have added a purely Gaussian 8T/T from the standard CDM model. The spherical 
imprint provided by the bubble is well evident in the bottom panel, but rather hidden in the top panel, 
depending of course on 8. Note however that both bubbles beco me strongly non-linear a t the present, and 
therefore are able to affect substantially the structure formation ( |Occhionero et al. 1997 ). The signal shows 



clearly its marked non-Gaussianity: the sound wave propagating outward and the central spot generated by 
SW and adiabatic effect determine the appearance of a few fascinating concentric rings of alternate colour 
on the sub-degree angular scale. As the mean anisotropy amplitude in Fig. (5) is not different from the 
ordinary Gaussian case, the distinctive imprints of these structures on the CMB are to sought not as much 
in the anisotropy amplitude ([l^), but rather in the phase correlation encoded in higher order statistics. 

The above results have to be extended to a realistic distribution of bubbles. The resulting maps will 
be compared with the present and next high resolution CMB o bservations. We leave these analysis to 
forthcoming works, beginning with the impact on the C; (|Io|), in Amendola et al. 1997. 



5. Conclusions 

The occurrence of first order phase transitions in the early universe is largely predicted, and its indirect 
detection through the present observable traces would be of crucial importance for our understanding of 
fundamental physics. If one or more of such processes occurs well before the end of inflation, the nucleated 
bubbles are stretched to large scales; this has very important consequences on the observability of these very 
high energy physics phenomena. Here we predict the traces left by inflationary bubbles on the CMB. To 
achieve this aim, we first model analytically a bubble. We use the general linear theory to get the Fourier 
transform of the corresponding perturbation in the photon baryon plasma and to evolve it until decoupling. 
Then, by considering a CMB photon decoupled at some point inside the bubble, we antitransform to get the 
8T/T carried by the photon; by repeating the above computation for each direction and for each decoupling 
point weighted with the last scattering probability function, we get the whole CMB anisotropy produced 
by an inflationary bubble. 

We find that the CMB perturbation propagates on the scale of the sound horizon at decoupling (~ 1° 
in the sky), and appears as a fascinating series of concentric rings of alternate color (different sign of 
8T/T), a general feature that remains true as the bubble is moved back and forth with respect to the last 
scattering surface. Therefore, in comparison with ordinary anisotropies from scale invariant perturbation 
spectra, the present signal presents a marked non-Gaussianity. The overall mean amplitude of the signal 
from a bubble with radius R and central density contrast 8 <C 1 is of the order of 8(R/ H^ 1 ) 2 , as expected 
for linear perturbations. Therefore, bubbles with comoving radii corresponding to the size of the observed 
large scale structures (tens of Mpc), and characterized by a central density contrast 8 ~ 10~ 3 , induce a level 
of anisotropy to be interestingly compared with the observed ones. 

If bubbles are produced during the inflationary era, the resulting fluctuation spectrum is a superposition 
of the standard Gaussian fluctuations from the slow rolling inflaton and the scale dependent non-Gaussian 
ones from the bubbles. The results of the present paper will be used to make quantitative predictions of the 
imprint of these inflationary models on the CMB. Since the angular scale involved by the bubbles are well 
below 1°, our computations shall be compared with the high resolution CMB maps provided by the next 
MAP and Planck experiments. 
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properties of the k-spectrum are almost the same in any of the case shown. Examples of symmetric shells 
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(cj = a a = a/2) arc on the left: a = .2 (short dashed line), a = .3 (solid line) and a = .5 (long dashed line) 
Cases of asymmetric shells are on the right: &i = .2, a D = .1 (short dashed line), <7j = <r = .15 (solid line) 
and cr, = .1, <7 = .2 (long dashed line). 

Fig. 2. — CMB perturbations at decoupling induced by an inflationary bubble: adiabatic effect in the top, 
Doppler in the middle, and Sachs- Wolfe in the bottom panel. The different lines refer to different radii as 
indicated. The Sachs- Wolfe effect does not exceed the bubble size, while the other effects generally extend 
on the sound horizon for the photon baryon plasma at decoupling. 

Fig. 3. — The outward traveling waves in the adiabatic (top) and Doppler (bottom) effects for a bubble with 
R = 20h~ 1 Mpc at different rcdshifts. The reduction of the amplitude is due to the Silk damping, that is 
becaming more and more active with time. 

Fig. 4. — CMB anisotropies induced by bubbles of interesting size near the LSS in different positions: d = —R 
(short dashed line), d = (solid line) and d = +R (long dashed line). The signal scales linearly with the 
central density contrast and involves angular scales up to 1°. 

Fig. 5. — A portion of microwave sky containing a purely Gaussian CDM ST/T field together with the 
signature of an inflationary bubble with R = 20/i _1 Mpc, 5 — .005, .01 (respectively top and bottom) and 
sitting exactly on the LSS (d — 0). Horizontal scale is in acminutes. Although the spherical imprint is well 
evident in the bottom panel but only hidden in the top one, both bubbles are non-linear by the present and 
capable to affect substantially the large scale structure. 
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